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A huge systematics of femtoscopic measurements have been used over the past 20 years to characterize the 
system created in heavy ion collisions. These measurements cover two orders of magnitude in energy, and with 
LHC beams imminent, this range will be extended by more than another order of magnitude. Here, I discuss 
theoretical expectations of femtoscopy of A+A and p + p collisions at the LHC, based on Boltzmann and 
hydrodynamic calculations, as well as on naive extrapolation of existing systematics. 
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I. INTRODUCTION 

What distinguishes ultrarelativistic heavy ion physics from 
particle physics is its focus on geometrically large systems. 
The desire is not to understand fundamental processes, as in 
the latter field, but to create and probe new states of matter 
and access the only phase transition associated with a funda- 
mental interaction (QCD). The geometrically-sensitive, bulk 
properties are the crucial ones, and these are reflected in the 
soft (p T ~ Aqcd) sector. 

In soft sector observables, long-term baselines have been 
established over a large energy range. Prior to first data at 
RHIC, it was commonly speculated (and hoped) that large de- 
viations from these systematics (e.g. tjP/k^ ratios, sidewards 
flow, strangeness enhancement, total multiplicity) would sig- 
nal clearly the qualitatively different nature of the system cre- 
ated there JU. In femtoscopic systems, rather generic ar- 
guments led to expectations J2l |3|] of a rapid increase, with 
y/siw, in the pion "HBT radii" R out and Ri on g, reflecting rela- 
tively long timescales of the transition from deconfined QGP 
to confined hadronic matter. 

Such dramatic speculations are largely absent today, in an- 
ticipation of LHC collisions. Soft-sector, global observables 
at RHIC are only quantitatively different than they are at lower 
energies. Even in the high-pr sector, where jet suppression 
and partonic energy loss measurements have generated huge 
excitement, energy scans at RHIC reveal that the data indicate 
more of an evolution than a revolution. 

This is all to the good. Discoveries via sharp jumps a la su- 
perconductivity are not our lot. The real science behind heavy 
ion measurements (at very high energies as well as at much 
lower ones) lies in understanding the details of the data. 

Femtoscopy 0], the geometric measurement of systems at 
the fermi scale, has been used to characterize the freezeout 
substructure of heavy ion collisions for two decades in time, 
and over two decades in collision energy [5]. Soon, this en- 
ergy range will be extended by another decade at the LHC. 

It will be important to understand the evolution of the non- 
trivial space-time substructure of the bulk system as the initial 
conditions change dramatically with energy. In Figure Q] is 
shown the number of refereed-journal papers of experimen- 
tal femtoscopic results in relativistic heavy ion collisions, as a 
function of publication year. Even armed with nothing more 
than this Figure and Reinhard Stock's observation |6[] that 
"HBT experiences a renaissance of new insights roughly ev- 
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FIG. 1: The number of refereed-journal publications reporting new 
femtoscopy results in relativistic heavy ion collisions. Beginning of 
availability of new beams are indicated in yellow boxes. 



ery five years," we may confidently expect a barrage of new 
femtoscopic information to digest within the next 5 years. 

Here, we discuss predictions for femtoscopy at the LHC. 
In the next Section, we consider the case of simple extrap- 
olation of measured femtoscopic trends, with no reference 
to physics per se. In Section [TTTJ we consider predictions of 
Boltzmann/cascade transport calculations, and in Section ITVl 
those of hydrodynamical models. Finally, we discuss spec- 
ulations on the physics behind femtoscopic measurements in 
p + p collisions, which will, in fact, be the first results avail- 
able at the LHC. At the end we summarize. 

II. NOTHING NEW UNDER THE SUN (NNUS) SCENARIO 

Femtoscopic measurements display rich, multidimensional 
and nontrivial systematic dependences upon kinematic (pr, 
y, etc) variables and particle species ||5l |7[]. The depen- 
dence upon global variables such as ^/snn and impact param- 
eter, however, appears significantly more trivial. Schemati- 
cally characterizing the measured femtoscopic length scales 
as a multidimensional function, evidence thus far indicates an 
overall factorization 
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FIG. 2: Pion HBT radii plotted versus the number of participating 
nucleons (left) and versus the charged particle multiplicity (right) for 
collisions of varying centrality and a wide range of energies. Com- 
pilation from @]. 



where is a dimensionless function containing, e.g. decreas- 
ing "HBT radii" with particle mj. The dimensional scale R g is 
determined by global observables. However, as indicated by 
the second equality of EquationQ] to good approximation the 
only relevant global observable is the total multiplicity M of 
the collision. In fact, this multiplicity dominance well apply 
to all soft-sector observables J^]. 

There are at least three caveats to the above statement. 
Firstly, the CERES [9] collaboration has shown that the scale 
depends also on the freeze-out chemistry (baryon-to-meson 
ratio) in addition to the multiplicity. This is important at 
low (AGS) energies. For collisions above top SPS energy, 
y/sNN ~ 17 GeV, the chemical evolution is sufficiently weak 
that one may consider multiplicity only. Secondly, as seen in 
Figure |2j there is some residual dependence of the outward 
radius on ■Jsnn in addition to multiplicity; in this Section, 
we ignore this potentially important detail. Finally, the az- 
imuthal (<|v, as determined relative to the reaction plane) de- 
pendence Hoi [Til [T2I [I3I1 likely at some point violates the 
factorization. While it has not been experimentally tested, 
two collisions producing the same multiplicity, one very pe- 
ripheral (i.e. spatially anisotropic in the entrance channel) 
at high energy and the other very central at low energy, pre- 
sumably generate freezeout distributions with different spatial 
anisotropy, which is then reflected in the azimuthally-sensitive 
femtoscopy lfl4ll . See SectionllVlfor further discussion. 

These caveats stated, however, the factorization of Equa- 
tion Q~|is probably our best, zero-new -physics guide to simple 
extrapolation of femtoscopic trends measured over two orders 
of magnitude in ^/saw and from from the lightest (p + p) to 
the heaviest {Pb + Pb) systems. Figure [2] suggests a simple 
form /? ? (M)ocM 1 /3 ; this ignores the finite offset ~ 1 fm when 



extrapolating M — ► 0, but this is negligible for high multiplic- 
ity. This relation may reflect that a constant freezeout density 
drives the femtoscopic scales ||9j], though this neglects any dy- 
namic effects. Assuming that this simple proportionality con- 
tinues, then, we know R g (M) and determining femtoscopic 
expectations boils down to anticipating the multiplicity at the 
LHC. 

A naive extrapolation (0. Ht] of systematics suggests that 
dN/dy at the LHC will be 60% larger than that observed at 
RHIC. Thus, the zeroth-order expectation is that length scales 
at the LHC will be 17% (1.6 l ' 3 = 1.17) larger than those 
measured at RHIC, for all kinematic selections and particle 
species, according to Equation[T] 

Going beyond simple extrapolation to include a physical 
picture, saturation-based calculations lfl5ll give much higher 
multiplicity- roughly triple that at RHIC. This leads to expec- 
tations of length scales 45% higher than those at RHIC. Thus, 
Riong for pions at midrapidity and low pj in central collisions 
would be 1 .45 x 7 fm = 10 fm. 

Multiplicity predictions based on Boltzmann/cascade cal- 
culations can be significantly higher yet. Selecting two 
for which femtoscopic predictions also exist (Section ITTTb , 
A Multi-Phase Transport (AMPT) calculation Jl^] and the 
Hadronic Rescattering Model (HRM) |T3l predict 5 x and 7x 
RHIC multiplicity, respectively. Thus, femtoscopic scales at 
LHC may be as much as 90% higher than at RHIC. Depending 
on the final-state interaction which produces the two-particle 
correlation function, measuring length scales of ~ 15 fm may 
challenge experimental two-track resolutions. For two-pion 
correlations, such scales are within the capabilities of the AL- 
ICE detector |fl8ll. 

HI. BOLTZMANN TRANSPORT CALCULATIONS 

More interesting than simple scaling relations are mod- 
els with real physics and dynamics, such as transport calcu- 
lations. Boltzmann/cascade transport models generally re- 
produce "HBT radii" at RHIC better than do hydrodynamic 
calculations [5]. The reasons behind this include different 
physics in the models, a more detailed description of the 
kinetic freezeout, and the use of more appropriate methods 
of calculating the radii [19]. Predictions of pion HBT radii 
with each of the transport calculations discussed in Section llll 
reveal predictions more subtle than the simple multiplicity- 
scaling discussed above. 

For an infinite and boost-invariant system (only an approx- 
imation of reality, of course), the longitudinal HBT radius 
Riong is proportional to the system evolution time (i.e. be- 
tween interpenetration of the ions and kinematic freezeout of 
the products) iHH^]. Naturally, this is not a unique, system- 
wide time, but a distribution. An example is seen in Figure [3] 
in which the pion freeze-out time distribution for collisions at 
RHIC and LHC are compared in the HRM calculation. The 
LHC timescales are roughly double those at RHIC. Although 
HRM is not explicitly a boost-invariant model, we see in Fig- 
ure |4] that Riong reflects this timescale increase, roughly dou- 
bling when the energy is increased from RHIC to LHC en- 
ergies. The ~ 70% increase in Ri ong is roughly consistent, 
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FIG. 3: The freezeout time distribution from the Hadronic Rescat- 
tering Model of Humanic fl7ll for RHIC and LHC conditions. 



then, with expectations from both a timescale and from the 
multiplicity-scaling point of view. This is certainly not a co- 
incidence, as the increased timescale is due in large part to the 
increased multiplicity. 

On the other hand, there is more going on. The predicted 
Pt -dependence of both Ri om , and are steeper at the LHC 
than at RHIC. Also, the increase in R S ide is significantly less 
than 90%. Both of these effects are consistent with a freezeout 
scenario with significantly increased transverse flow IU411 . In- 
deed, transverse momentum distributions predicted by HRM 
are significantly harder (less steep) at the LHC than those at 
RHIC. Since the pr dependence of HBT radii |5j] and spectra 
in the soft sector are observed to change very little between 
y?SNN =20-i- 200 Ge V, it will be interesting to see whether 
this trend is broken at the LHC, as predicted by HRM. 

The HRM model is a deliberate effort to use the simplest 
(often criticized as too simplistic) physics picture, free of 
novel phases like QGR It is a pure hadron-based transport 
calculation, though the initial conditions may be taken from 
Pythia or Saturation-based scenarios 11711 . On the other side 
of the "simplicity spectrum" is AMPT, an attempt to describe 
the various stages of the system's evolution in terms of the 
most appropriate model for that stage Ir2lll . 

Similar to HRM, AMPT predicts stronger transverse flow 
at the LHC, as compared to RHIC, leading to steeper pj- 
dependence of HBT radii. In terms of scale, the transverse 
(longitudinal) radii are predicted to increase by 10% (30%). 
This is more modest than the predictions of HRM (30% 
and 70%, respectively), and much more modest than pure- 
multiplicity scalings of Sectionllll 

Thus the dynamical physics, in these models, lead to ex- 
pected details significantly beyond simple extrapolation of 
lower-energy results. 

IV. HYDRODYNAMICAL CALCULATIONS 

As mentioned, hydrodynamical models tend to repro- 
duce femtoscopic measurements more poorly than do Boltz- 
mann/cascade calculations. On the other hand, they have en- 
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FIG. 4: HBT radii from fits to pion correlation functions from the 
Hadronic Rescattering Model of Humanic 11711 for RHIC and LHC 
conditions. 



joyed huge success in reproducing momentum-space observ- 
ables such as elliptic flow. Furthermore, the conditions at 
LHC are likely to provide an even better approximation than 
at RHIC to the zero-mean-free-path assumptions of pure hy- 
drodynamics. Finally, the direct connection between hydro- 
dynamics and the Equation of State of strongly-interacting 
matter (color-confined or not) remains a compelling reason 
to explore soft-sector, bulk consequences of the model. 

A. Source Length Scales 

Recently, Eskola and collaborators l22tl coupled a 
pQCD+saturation-based prediction for initial conditions at 
LHC to their 1 + 1 -dimensional hydro calculation. The 
Equation-of-State featured a first-order phase transition be- 
tween an ideal QGP at high temperature and a hadron reso- 
nance gas at low temperature. 

As shown in Figure[3] the initial energy density at which hy- 
drodynamics is assumed to take over expected to roughly an 
order of magnitude larger at the LHC than at RHIC, due both 
to increased gluon production and to shorter system formation 
(thermalization) time To at the higher energy. Since the initial 
transverse scale changes only little, the pressure gradients will 
likewise be much higher at LHC, leading to increased trans- 
verse flow. These effects place competing pressures on the 
space-time evolution of the system, and on the femtoscopic 
scales at freezeout, as discussed below. 

The increased energy density (directly associated with en- 
tropy density and thus multiplicity) tends to produce longer 
timescales at the LHC. Longitudinal expansion tends to cool 
the system towards freezeout conditions. However, especially 
at the LHC, the large transverse flow generated by the intense 
pressure gradients cannot be ignored. Eskola and collabora- 
tors 112211 . estimate that the time required to cool from the max- 
imum energy density (at r = in Figure [5] to the critical en- 
ergy density (e e = \.92>GeV / fm i ) would be 6 fm/c (20 fm/c) 
at RHIC (LHC), due to longitudinal expansion alone. How- 
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FIG. 5: Initial energy density distribution in the transverse plane, 
calculated by Eskola et al l22ll . 



ever, when transverse dynamics are included, the cooling 
times become 5 fm/c (7.5 fm/c) at RHIC (LHC). The evo- 
lution time to kinematic freezeout- say until T w 140 MeV- 
is To ~ 12 — 14 fm/c in both cases; this is the timescale most 
directly probed by femtoscopy. This is dramatic- the effect 
of transverse flow on cooling timescales can almost be ne- 
glected at RHIC, while it is dominant at the LHC. This is rem- 
iniscent of the cascade calculations discussed in Section Hill 
the much stronger flow may well lead to deviations from the 
trends (e.g. pj -dependence of pion HBT radii being indepen- 
dent of ^/snn) established so far at lower energy. This aspect 
of the NNUS scenario may finally be violated. The qualita- 
tive difference is apparent from the freezeout hypersurfaces 
at RHIC and LHC, shown in Figure [6] The Figure is from a 
calculation by Kolb and Heinz |23y, but is similar to Eskola's. 

It would be very interesting to know whether the other as- 
pect of NNUS, namely the multiplicity scaling shown in Fig- 
ure |2] is satisfied by the hydro models. Unfortunately, Es- 
kola did not calculate pion "HBT radii," and Heinz and Kolb 
did not calculate multiplicity, so a consistent estimate of the 
scaling cannot be checked. However, the former predict that 
the multiplicity at LHC will be approximately triple that at 
RHIC, corresponding to a 40% increase in HBT radius un- 
der NNUS scaling. Heinz and Kolb B23I1 do, in fact, predict 
roughly this increase in the transverse radii, but- remarkably- 
they also predict a significant decrease in Ri ong at LHC rela- 
tive to RHIC! S 



B. Source Shape 

Azimuthally-sensitive pion interferometry- the measure- 
ment of spatial scales as a function of emission angle relative 
to the reaction plane- probes the shape of the freezeout con- 
figuration, in addition to its scale lfTollT4ll24ll . At finite impact 
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FIG. 6: Left: freeze-out space-time hypersurface in collisions with 
finite impact parameter calculated by Heinz and Kolb 12311 . The 
"IPES" calculation is an estimate of the system created at the LHC. 
The impact parameter is defined to lie in the x-direction. Right: 
time-,z- and momentum-integrated freeze-out shapes in the trans- 
verse plane. 



parameter, both the spatial configuration of the entrance chan- 
nel and the resulting momentum distribution in the exit chan- 
nel are anisotropic. In particular, the initial state is spatially 
extended out of the reaction plane, and the resulting flow is 
stronger in the reaction plane (elliptic flow v>2 > 0). 

Due to the preferential in-plane expansion, as the system 
evolves the spatial configuration should become increasingly 
in-plane extended (equivalently, decreasingly out-of-plane ex- 
tended). Thus, knowledge of the entrance-channel shape 
(e.g. though Glauber model calculations) and measurement of 
the exit-channel shape (through femtoscopy) provide "bound- 
ary conditions" on the dynamical spacetime evolution of the 
anisotropic system, and probe the evolution timescale. The 
extracted timescale is model-dependent, requiring in princi- 
ple a detailed time evolution of the flow. However, a simple 
estimate B25I1 of the timescale extracted through shape mea- 
surements and that extracted from blast-wave fits 13, |H to 
azimuthally-integrated HBT radii are roughly consistent. 

Measurements of the freezeout shape at the AGS Hill and 
RHIC 1IT2I [I3I1 indicate an out-of-plane-extended configura- 
tion. Consistent with the fact that preferential in-plane expan- 
sion (i.e. elliptic flow) is stronger at RHIC, the configuration 
at the higher bombarding energy is rounder. This is shown in 
Figure |7J in which the transverse anisotropy is characterized 
by £ = (R? - R\)/(Ry + R%) (x is in the reaction plane). 

The anisotropic shapes and corresponding azimuthally- 
selected HBT radii have been calculated in two dynamical 
models. At the AGS, the transport code RQMD||Z7| repro- 
duces the overall scale H28I1 and the anisotropy fl 1 0L ll 111 of the 
source reasonably well, as shown in Figure [7] At RHIC, the 
2 + 1 hydro code i23ll reproduces the shape quite well, while 
missing the scale. At the LHC, the latter calculation predicts- 
again- a qualitative change in the freezeout distribution. As 
shown in the right panel of Figure [6] the source is expected to 
evolve to an in-plane configuration (e < 0). 

However, again, this huge flow has other dramatic and qual- 
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itatively new implications: the decreased Ri on g mentioned in 
Section llV Al and an actual sign inversion of the oscillations 
of the HBT radii with p T ITIB . 
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much weaker pj -dependence of Rj nv . 



V. PROTON COLLISIONS 

Before heavy ions are accelerated at the LHC, proton col- 
lisions at y/s = 1.4 TeV will be measured. While the thrust 
of the p + p program is towards Higgs physics, it is well- 
recognized that p + p collisions serve as a valuable reference 
to heavy ion analyses in the "hard" (high-p^) sector, where 
one looks for the effects of the medium on particles coming 
from well-calibrated fundamental processes. 

Soft-sector analyses, too, should be performed for systems 
from the smallest to the largest, and the results compared. 
Since such analyses are assumed to measure the bulk prop- 
erties, one might well hope for qualitative differences when 
comparing results for p + p to Pb + Pb collisions. 

While pion HBT measurements have been common in both 
the high-energy and heavy-ion communities for many years, 
a direct "apples-to-apples" comparison between results from 
A+A and p + p collisions has not been possible until very 
recently. The STAR Collaboration at RHIC has reported the 
first direct comparison of pion HBT radii in Au+Au, Cu+Cu, 
d+Au, and p+p collisions, using the same detector, same en- 
ergy, identical techniques (event mixing, etc) to create the 
correlation function, identical coordinate systems and iden- 
tical fitting techniques 02911 . Remarkably, Gaussian fits to 
the correlation functions return "HBT radii" which factor- 
ize according to Equation [T] i.e. Fj, which quantifies the 
dynamically-generated substructure, is identical in p + p and 
A+A collisions. However, the STAR data show significant 
non-femtoscopic structures B29I1 . which must be properly ac- 
counted for [30] before drawing firm conclusions. If the fac- 
torization is unchanged after a more sophisticated treatment, 



the physics implications might be dramatic. We do not dis- 
cuss this here, but simply observe that the NNUS scenario is 
a likely baseline expectation for p + p at LHC. 

The increase of HBT radii with multiplicity has also been 
observed previously in p—p collisions by the E735 Collabo- 
ration [32]. While in A +A collisions, this is naturally related 
to increasing length scales in the entrance channel geometry, 
Paid and Skowrohski [33] postulate that jet dynamics, rather 
than bulk properties, drive this dependence in the p + p sys- 
tem. Within a simple model of hadronization, they can repro- 
duce the E735 multiplicity dependence, and make predictions 
of similar multiplicity dependence for p + p collisions at the 
top LHC energy. However, since the expectation of increas- 
ing length scales with multiplicity seems to be rather generic 
to all scenarios, it will be interesting to see these predictions 
expanded to more differential measures- say, the multiplic- 
ity and pj dependence, probing both aspects of Equation Q] 
This should allow a more discriminating comparison between 
models, allowing some to be ruled out. 

Such differential predictions have very recently 
been performed by Humanic Ir3lll . in the context of a 
Pythia+hadronic rescattering (through HRM) scenario. It is 
found that, contrary to some expectations, hadronic rescat- 
tering is crucial to understand the M ® pj dependence of the 
E735 data. Reproducing the data requires the assumption 
of a surprisingly short hadronization timescale (~ 0.1 fm); 
longer timescales do not allow sufficient hadronic rescattering 
needed to describe the /^-dependence. The prediction of 
the model for the highest multiplicity p + p collisions at 
xfs = 1 .4 TeV are shown in Figure [8] 
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VI. SUMMARY 

Two decades' worth of femtoscopic systematics Jft] in 
heavy ion collisions reveals a strikingly consistent and simple 
structure. The kinematic and particle-species dependences, 
which reflect dynamic substructure, decouple from the global 
scale, which depends (almost) solely on multiplicity. The as- 
sumption that the factorization of Equation[T]persists- i.e. that 
Fk remains unchanged at the LHC- together with the assump- 
tion R g {M) ~ \/M, is the essence of the NNUS scenario. 

In the NNUS picture, all femtoscopic length scales mea- 
sured at RHIC will be reproduced at LHC, only scaled up 
by 20-90%, depending on the multiplicity prediction. Thus, 
the "pion HBT radius" Rt ong at low pj might be expected in 
the range (1.2 -=-1.9) X (7fm) = 8.4 -f- 13fm, while the aver- 
age shift between pions and kaons (about 6 fm at RHIC (H3l) 
would be in the range 7.2 -f- 11.5 fm. However, dynamical 
models generally predict interesting violations of NNUS at 
the LHC. The significant dispersion between predictions holds 
out the possibility that the data will eliminate some models. 

In the HRM and AMPT Boltzmann/cascade calculations, 
increased rescattering due to the higher density at the LHC 
generates much stronger global space-momentum correla- 
tions. This leads to flatter pj spectra for high mass particles, 
and to a steeper /^-dependence of the femtoscopic length 
scales; that is, F*. would pick up a yfs dependence, violating 
NNUS factorization. Transverse (longitudinal) scales are ex- 
pected to increase 10-30% (30-70%), relative to RHIC values. 

In hydrodyanamical calculations, much higher energy den- 
sities and pressure gradients at the LHC may generate quali- 
tatively new femtoscopic signals. Contrary to the situation at 
RHIC, the transversely explosive nature of the source at the 
LHC severely shortens the time until freezeout. The freezeout 
hypersurface is of a qualitatively different shape in transverse 
position and time; while transverse radii may increase by 40% 



relative to RHIC values, the longitudinal ones should expand 
little, and may even decrease. 

The evolution timescale may also be probed by measuring 
the anisotropic shape of the source in coordinate space, for 
non-central collisions. Here again, a qualitative difference is 
predicted by hydrodynamics between RHIC and LHC colli- 
sions. In particular, the greatly increased flow and somewhat 
increased evolution time lead to predictions of an in-plane ex- 
tended source, producing HBT radius oscillations 180° out of 
phase with those seen at lower energy. 

Probably as important as soft-physics analyses in heavy ion 
systems are parallel ones for p + p collisions. First prelim- 
inary "apples-to-apples" comparisons of Gaussian HBT ra- 
dius measurements at RHIC suggest that NNUS factorization 
continues to hold even for these smallest systems; it remains 
to be seen whether this conclusion survives more sophisti- 
cated treatment of non-femtoscopic correlations in the data, 
presently underway. In the context of two simple models, pion 
HBT radii at the LHC depend strongly on the hadronization 
scenario. Both predict an increase in femtoscopic freezeout 
scales with increasing multiplicity, which in itself will not dis- 
tinguish these models from any other. However, in one, the pj 
dependence is found to depend strongly on the hadronization 
time and degree of subsequent hadronic scattering. Such scat- 
tering is usually ignored in treatments of p + p collisions; in- 
deed, the lack of significant rescattering is believed to be their 
primary virtue as a reference measurement. As hinted at by 
the first RHIC measurements, maybe they are not so different 
from A+A collisions after all. More detailed measurements 
at the LHC may, in fact, spur a re-evaluation of ideas of the 
spacetime evolution of both heavy ion and hadronic collisions. 

In any case, we may confidently expect considerable activ- 
ity and excitement as the next mountain forms on Figure Q] 
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